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Abstract 

We study the matching of a general spherically symmetric spacetime with 
a Vaidya-Reissner-Nordstrom solution. To that end, we study the proper- 
ties of spherically symmetric electromagnetic fields and develop the proper 
gravitational and electromagnetic junction conditions. We prove that generic 
spacetimes can be matched to a Vaidya-Reissner-Nordstrom solution or one 
of its specializations, and that these matchings have clear physical interpre- 
tations. Furthermore, the non-spacelike nature of the matching hypersurface 
is proved under very general hypotheses. We obtain the fundamental result 
that any spherically symmetric body, be it in evolution or not, has un upper 
limit for the total net electric charge that carries. 
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1 Introduction 

Frequently, within the framework of General Relativity, in order to describe new 
or interesting physical situations, or to investigate certain theoretical aspects, the 
construction of a model is required. This is commonly done by matching two given 
and known space-times. Thus, for instance, models for stars (to study the genesis of 
black holes, the possibility of escaping a black hole region, the violation of the cosmic 
censorship conjecture, etc.) and models for local inhomogeneities in a cosmological 
context (to study the origin of galaxies in the primitive Universe, the possibility of 
primordial black holes, the evolution of the stars in the Cosmos, etc.) can be studied 
by means of the junction of two different space-times. 

With regard to the first type -construction of stellar models- the common fea- 
ture in most of them is the matching of two space-times with special characteristics: 
one of them, which describes the so called interior of the star; and another, with- 
out matter but possibly with electromagnetic fields and/or null radiation, which 
represents its exterior. Both regions are separated by a timelike hypersurface, rep- 
resenting the surface of the star. A very typical simplifying assumption is that 
the complete space-time possesses spherical symmetry, so that in order to describe 
the exterior of an object characterized only by its mass the Schwarzschild vacuum 
solution is used; if, in addition, the space-time possesses a net electrical charge, 
then the Reissner-Nordstrom (RN) solution is used; and if radial null radiation is 
included then, instead of the previously cited metrics, the solutions of Vaidya 
and Vaidya- Reissner-Nordstrom (V-RN) fl30| , are used, respectively. 



Thus, we can classify the works for stellar models according to their exterior. 
The archetype with a Schwarzschild exterior is the pioneering work by Oppenheimer- 
Snyder [f£|, in which the collapse of a dust cloud was studied. This work was later 
generalized from two different points of view: on the one hand, by considering a 
more general interior matter distribution. This was done for instance by Misner and 
Sharp p9| , where a perfect fluid interior was considered, and by Bel and Hamoui 



for the case of a fluid with anisotropic pressures. On the other hand, by allowing 
for the existence of an electromagnetic field. For example, when a charged dust 
interior with a Reissner-Nordstrom exterior was considered. Outstanding examples 
of this treatment were the works by Novikov, de la Cruz & Israel, Markov & Frolov, 



Hamoui, Vickers, Misra & Srivastava, Raychaudhuri and de Felice Sz, Maeda HI 
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Recently, Amos Ori [44| has found the explicit general solution for the Einstein- 
Maxwell's equations in this case. 

The next natural step for stellar models with a RN exterior was to consider 
an even more general interior. This was initially achieved by Bekenstein fl|], who 
considered a charged perfect fluid interior, by using a generalization of the formalism 
(for neutral perfect fluids) due to Misner and Sharp. Some other interiors were 
later considered. This is the case of the article by Florides |T!|, in which an exact 
solution of the Einstein-Maxwell equations for a static, spherically symmetric charge 
distribution is found. This represents an extension, to the charged case, of a previous 



solution by Synge |49|| . 

With regard to the radiating stars, that is, models with a Vaidya exterior, we have 
to point out the works by Herrera et al. (see for a summary of their results) and 
Santos et al. (see and references therein), where the interior energy- momentum 
tensor is interpreted as a fluid with heat flux. However, these works study only the 



shear-free fluid case. Less restrictive is the article by Fayos et al. [[16] where the 
preceding works on Schwarzschild and Vaidya [B2] interiors were generalized. 

The complete case with a radiating and charged exterior (that is, a V-RN exte- 



rior) has been treated by Oliveira and Santos ||42|| , but only in the case of a charged 
fluid interior with heat flux but without shear. This article represents a general- 
ization of the work by Misner and Sharp, and also, of the one by Bekenstein, in 
which dynamical equations for the spherical collapse are given. Likewise, Herrera 
and Nunez considered a V-RN exterior, but only for a charged shear- free perfect 
fluid interior, and with a flux of non-polarized radial radiation. In their article they 
extend the HJR method -described in this same reference- in order to obtain the 



evolution of radiating charged fluid spheres [ 3SJ . They use a heuristic ansatz which 
allows them to integrate the Einstein-Maxwell equations taking into account the 
matching conditions. 

Regarding the construction of local inhomogeneities in a cosmo logical context, 
the first fundamental work is the historically relevant paper by Einstein and Strauss 

, in which "the influence of the expansion of space on the gravitation fields sur- 
rounding the individual stars" was studied. They considered a vacuum Schwarzschild 
interior surrounded by a Universe, modelized by a Robert son- Walker exterior. The 
generalization, replacing the Schwarzschild solution with a Vaidya interior, in order 
to describe the so called primordial black holes, can be found in |21], [IT], [31], [15 . 
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Lake and Hellaby |H| used these models to study the formation of naked singu- 
larities in the collapse of radiating stars, as candidates for counterexamples of the 
cosmic censorship conjecture. 

It should be emphasized here that the actual matching procedure for the Einstein- 
Straus and Oppenheimer-Snyder models are mathematically identical and indis- 
tinguishable, because the two glued spacetimes are exactly the same: vacuum 
Schwarzschild with dust Robertson- Walker. They are in fact an outstanding ex- 
ample of complementary matchings, see |1| . This is in fact a general feature of the 
matching procedure and thus we can always attack the two problems (description of 
stars, and modelization of voids) in a common and unique framework. This is one 
of our aims in this paper, and to that end we will follow the philosophy put forward 
18| which will allow us to treat the two possibilities jointly in a natural way. 



in 



As a second goal of our paper, we wish to treat the fully general case in which 
the matter content in the interior for stars (or the exterior for voids) is left unre- 
stricted, and the corresponding exterior (resp. interior) can contain both radiation 
and electromagnetic fields. We will attack this problem by studying the matching 
of two spherically symmetric space-times V and V through a time-like matching 
hypersurface £ in such a way that one of the space-times will be completely general 
(V), and the other (V) may contain an electromagnetic field and/or radial null radi- 
ation. This means that V will be described by a V-RN solution or one of its special 
cases. We remark that our results will be general and independent of the different 
combinations that might exist in V: vacuum, only electromagnetic field, only radial 
null radiation, or radiation and electromagnetic field. 

The possible physical justification for including electromagnetic charge in the 
models arises from several different considerations: i) first of all, in the formation of a 
stellar object the repulsive Coulomb force acting in every charged particle of the same 
sign of the net charge of the star cannot be arbitrarily big, as is obvious. A simple 
calculation, taking into account the opposing forces acting inside a star with 
baryons, indicates that the net charge Q must obey the inequality Q < lCT 36 ^ f20|| . 
This means that only a small value of the net charge per nucleus is permitted, but not 
necessarily zero; ii) the existence of an electromagnetic field has drastic consequences 
on the global space-time structure -compare, for instance, the Schwarzschild solution 
with the Reissner- Nordstrom solution-. Thus, the evolution of a star with a net 
electric charge Q could be totally different to the evolution of a neutral star, no 
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matter how small Q may be. New possibilities -like that of a star avoiding the 
collapse towards a singularity- have been investigated since the appearance of the 
previously cited article by de la Cruz and Israel, in which it was shown for the first 
time that a charged object, after collapsing beyond the event horizon, can expand 
and re-emerge in another asymptotically flat region. Of course, all this depends on 



the stability of the Cauchy horizons, see e.g |50| and references therein, a problem 
which is still subject to controversy; finally, even if there is no electromagnetic field 
in V, we could allow for a radial charge distribution in V with total net charge 
vanishing. This is necessary to describe charge redistributions that may take place 
in the interior of a star. A phenomenon like this happens, for instance, in the hybrid 
stars, in which there is global, but not local, charge neutrality, see |20 . 



The plan of the paper is as follows. First, we will study, in section ^|, the 
properties of general spherically symmetric space-times. In particular we write down 
the necessary and sufficient conditions which guarantee the absence of curvature 
singularities. Since we want our models to have a physical meaning, we impose the 
fulfillment of the dominant energy conditions and, as a first consequence of this, we 
obtain the conditions such the mass function will be non-negative. In section |3| we 
analyze the electromagnetic fields, both null and non-null, that are compatible with 
the spherical symmetry, and the restrictions on the distribution of charge that arise 
at this level. Then we devote section ^ to the study of the spherically symmetric 
V-RN space-time V. In section [|, the gravitational and electromagnetic matching 
conditions are derived. We distinguish the cases with or without a surface charge 
current density in the matching hypersurface. At the end of this section we study 
the matching results, their interpretations and some consequences. Next, in section 
|6], we analyze the restrictions on the matter contents of the space-time V derived 
from the matching conditions. In order to guarantee the physical meaning of V 
-no matter which interpretation of the matter content is used-, in section [7| we will 
demand the fulfillment of the dominant energy conditions to the energy-momentum 
tensor, and we will examine the physical consequences that such constraints imply. 
We end up with some conclusions, in which the main results are summarized: they 
are the existence of a maximum for the total charge of the bodies, the matchability 
of generic spacetimes to V-RN or its specializations -these matchings having a clear 
physical interpretation — , and the causality of the matching hypersurface in general. 
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2 Basics on spherically symmetric space-times. 

Let us consider a four- dimensional spherically symmetric space-time V, so that its 
line-element can be expressed in radiative coordinates j|] = {u, R, 8, ip} (// = 
0,1,2,3) as 

ds 2 = -e AI3 xdu 2 + 2ee w du dR + R 2 rffi 2 , (1) 

where x = l — 2m/R, e 2 = 1, (3 and m depend on {u, R}, and dfl 2 = d9 2 + sin 2 9d<p 2 . 
It is easily checked that m(u,R) is the well-known mass function defined by Cahill & 
McVittie ||, which represents the total energy inside the two-spheres with constant 
values of u and R, see e.g. p3| . 

The global structure and general properties of the space-time depend on the 
behavior of the functions m(u, R) and /3(u, R). In this sense, for example, it can be 
seen that no spherical closed trapped surfaces (see |22[ for definitions) are present 
in the whole spacetime if and only if \ is non-negative everywhere (that is to say, 
2m(u,R) < R). In fact, the two-spheres defined by u = u c =const. and R = 
R c =const. are closed trapped surfaces if and only if x( u c R c ) < 0> an d thus the 
hypersurface x — 0, which is the boundary between regions with and without closed 
trapped spheres, is the apparent horizon. 

Similarly, all curvature invariants will be finite at R — |5T| if and only if 
(comma denotes partial derivative) 

# fl (0) = ; m(0) = m, R (0) = m, M (0) = 

so that we have 

Proposition 2.1 The necessary and sufficient conditions preventing the existence 
of a curvature singularity at R = are 

v P(u,R)-l3o(u) 
m(u,R) 

lim — = m 3 {u), (2) 
-R^O R 6 

where (3q{u) = lim^o (3(u, R), ^(u) and m^{u) are finite functions of u. 

We shall later use these conditions. 

In order to build up a model, some knowledge or constraints on the energy- 
momentum content, evaluated over the region of the space-time under study, are 
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compulsory. The standard energy conditions, see e.g. [22|, are among these physical 
constraints, so that in this paper we will always demand that the space-times under 
consideration fulfill the dominant energy condition (DEC), which requires that the 
energy-flux 4- vector is always non-spacelike ||22|| . For spherically symmetric space- 
times with the metric expressed in radiative coordinates, the DEC reduces to the 
following inequalities: 

P,r > 0, (3) 

m, R -R X [3, R > 0, (4) 
2e" 2/3 

A > \P 2 \, (6) 



where 



1/2 



2\ 2 2 ( 2e~ 2/3 

3 X 6 1 

A = ^ + 4x(3,R-—m, R )(3,R+2x(3,RR-—m,RR+e2e- 2 P(3, U R. 

The inequality (f|) can be rewritten as 

(me 2 ^),R > ^), R > 

where the last inequality follows from (|3|). Take an arbitrary null hypersurface 
u = w c =const. By integrating the above inequality in such a hypersurface between 
Ri and R 2 (> R±) we obtain 

„R 2 D 

m{u C) R 2 )e 2 ^ u ^ - m{u ci R^e 2 ^^ > J r ^( e 2 ^ u ^) tR dR > 
so that 

m(u c , R 2 )e w{Uc ' R2) > m(u c , Rje 2 ^'^ (7) 

From here we deduce that if m(u c ,Ri) > 0, then m(u c ,R) > for all R > R\. In 
particular, it is sufficient that \imR_ m(u c , R)e 2/3( - Uc ' R ^ > for the mass function to 
be non-negative on the whole hypersurface u = u c . Taking into account that u = u c 
is arbitrary, we have proven the following 

Proposition 2.2 // a spherically symmetric space-time is such that the mass func- 
tion m(u, R) is non-negative at R = and the dominant energy conditions are 
fulfilled, then m(u, R) > for all R > 0. 



3 ELECTROMAGNETIC FIELDS WITH SPHERICAL SYMMETRY 8 



Actually, not all DEC's are needed here, but just conditions An interesting 

case for this proposition, taking into account fl2|), is the following (compare with the 
more restrictive prop. 6 in p3| ) 

Corollary 2.1 If a spherically symmetric space-time has no curvature singularity 
at R = and the DEC are fulfilled, then the mass function m{u, R) is non-negative 
everywhere. 

3 Electromagnetic fields with spherical symmetry 

Starting with a spherically symmetric space-time in its form (p]) we can compute, 
through Einstein's equations, its energy-momentum tensor T. Choosing an or- 
thonormal cobasis {V^} = {u, n, u;g, u) v } such that ujq = R d9, uj^ = R sinO dip, 
n is a space-like 1-form orthogonal to ujg and u> v , but otherwise arbitrary, and u is 
a time-like 1-form which completes the orthonormal tetrad, one gets: 

T = Au<g>u + 5n<g>n + C(u<g>n + n<g>u) + 

+D (u)0 <g> uj e + uj v <g> u v ) (8) 

where A, B, C and D are functions of u and R. This energy-momentum tensor will 
be, in general, a mixture of fluids, gas, radiation and perhaps an electromagnetic 
field. For the sake of generality, we shall not identify all the particular fields and 
components of the matter contents. Nevertheless, avoiding this general rule, we 
shall in fact identify the electromagnetic part of T. This will reveal very useful 
because, among other things, it provides junction conditions for the electromagnetic 
field itself, independent of the gravitational ones, a fact which plays a central role 
in the physics of the problem. Thus, throughout this paper we will speak of the 
spherically symmetric electromagnetic field F, which is of course a 2-form solution of 
the Maxwell equations and must be considered as a partial source of the gravitational 
field so that 

T = P + E (9) 

where E is the part associated to the electromagnetic field 

E = ^-{F^Ff - -F lS F^ Va p) V a ® V 3 (10) 

and P is the non- electromagnetic part. Obviously the case without electromagnetic 
field is trivially included here when F = 0. Assuming that P also has the structure 
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(H), then E will share this same structure, so that T, P and E will all be invariant 
under the action of the SO (3) group. 

Consider a general 2-form F and demand that its energy-momentum tensor flT0| ) 
has the structure (H). This imposes several restrictions on its components and one 
can easily prove that only two possibilities arise, depending on whether F is null or 
not: 

F = ijj u A n + (f) (jl> 9 A Up = & (cos 1] u A n + sin rj ujq A uj^) , (non-null) (11) 
F = $ [cos 77(1 A uj ) + sin 77 (1 A (null) (12) 

where $ = $(«, R) and r\ = r](u, R, 9, if), and 1 is a future-pointing radial null vector 
field tangent to the (outgoing or ingoing) radial null geodesies. We also assume that 
1 is affinely parametrized. The corresponding energy-momentum tensors E read 



$2 

8~7T 

$2 



E = — (u <g> u - n <g> n + uj q <g> u) + uj^ <g> u) 9 ), (non-null) (13) 



E = — 1 ® 1 . (null) (14) 

47T 

The previous restrictions are purely algebraic. In addition, the electromagnetic 
field F must also satisfy the Maxwell equations 

dF = (15) 
SF = 4ttJ (16) 

where 5 is the co-differential (or divergence operator) and J is the electromagnetic 
current 1-form. These equations impose further restrictions on the components of 
F as follows: 

a) Non-null electromagnetic field. Equations ([T~5| - [T6D applied to the electromag- 
netic field (HTD lead to 



d(R 2 (f)) = (j) = c/R 2 (17) 

_J [i(u)d(R 2 il))n + i(n)d(R 2 il>)u] = J (18) 



4ttR 2> 

where c is a constant[|, and z() indicates inner contraction. This implies that 77 (and 

therefore <fi and ip) can only depend on u and R. 

1 The possibility c ^ will be constrained later in this section and will turn out to be impossible 
in our models as proved in section fj. 
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Notice that J can be split into two components, one along u and the other along 
n. Thus, an observer whose world lines are tangent to u will measure a charge 
density given by 

i(n)d(R^) 
Q = - J Ua = An& 

and a conduction current || 

i(u)d(R^) 

J = J-£U = ^^^n. (19) 

The total charge q within a 2-sphere S defined by constant values of u and R is 0] : 

g=i-/ s *F =► g = ^ 2 (20) 

where *F is the Hodge dual of F. 

b) Null electromagnetic field: The Maxwell equations ( |I7|]T6| ) applied to the 
electromagnetic field flT2"|) provide now 

i(k)d(R<$>cosr]) = , (21) 
i(k)d(i?$sinr/) = , (22) 
i(u;^)d(cos77) sin^ — 2(<^)d(sinr/sin6>) = , (23) 

$ 

— [i(u}g)d(cosr] sin 0) esc 9 — i(uj v )d(smr])] 1 = J , (24) 
An 



where k ■ k = 0, 1 • k = — 1 and, in order to write J in the form (^J), we have used 
equations ( |2~T| - |2"2"D . From here we infer that the electromagnetic 4-current ^ can only 
be null or zero and that the charge in any 2-sphere must be zero: q = 0, V(ii, R). 
Observe that if smi] / we can write (0), using (|23|), as 

$ 

J = esc i] i(u}g)drj 1 (25) 

An 



while if sin rj = then ( p4|) becomes 

$ 

J = cot 6 1 . (26) 

The case of physical relevance is that with no sources for the null electromagnetic 
field or, in other words, with J = 0. From ( Po] - PBD , we see that this is only possible 
if i(ujg)dri = and sin 77 7^ 0. However, (|23|) implies then that r), v = cos8, which is 
clearly incompatible with the previous formulas. Therefore, we conclude that there 
cannot exist a sourceless null electromagnetic field whose energy-momentum tensor 
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is SO (3) invariant. If these sources existed the associated electromagnetic 4-current 
should be a null vector, a possibility which has a questionable physical interpretation 
(see e.g. ||, QTTJ and references therein). Thus, let us remark the well-known result 
that the only electromagnetic field invariant under the action of the group 50(3) is 
the non-null one defined by (|lli|l7| - ^) . 

As is well known, the tensor (|13"|) obeys the dominant energy conditions by 
itself, and we have assumed that the whole energy-momentum tensor (||) so does. 
Nevertheless, this does not imply that P must satisfy the DEC. However, as we 
wish to avoid negative energy densities or space-like energy fluxes associated with 
P, in what follows we are going to demand the fulfillment of the dominant energy 



conditions for this part too. Taking into account (0), a straightforward calculation 
proves that conditions (|^) and @ remain the same while the inequalities (|]) and 
(Q) must be replaced by the stronger requirements: 

$ 2 < -^(m, R -RxP, R ), (27) 
2$ 2 < A + P 2 , (28) 
< A-P 2 . (29) 
At this stage, we can already deduce important consequences from these inequal- 



ities. For instance, the energy condition (]27 ) has an important physical interpreta- 
tion. Using that <3> 2 = 4> 2 + ip 2 , the definition ( |2"0"| ) of q, and the result ( |TTD for 0, we 
can rewrite (|27D as: 

q 2 + c 2 <2R 2 (m, R -R X (3, R ). (30) 

Recall that q is a function of u and R. Thus, if the space-time is given (that is, if 
we know m(u,R) and f3(u,R) explicitly), then this relation can be interpreted as 
providing a maximum value to the charge enclosed in every 2-sphere. 

If we now demand also the absence of curvature singularities at R = (see (Q)), 
from (^) we get that, in a neighborhood of R = 0, 

q 2 + c 2 < 2[3m 3 (u) — 2j3 2 (u)]R 4 + higher order terms. (31) 

Hence, as c is a constant and the second member of the inequality is zero for R = 0, 
we get that c = and q(u, 0) = 0, and by using (|17D , we arrive at the following 

Proposition 3.1 If a spherically symmetric space-time has no curvature singular- 



ities at R = and the energy condition (21) is fulfilled, then 

= 0, q(u,R^ 0) = 0. (32) 
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In other words, the radial observers will not measure any magnetic field (since its 
value is proportional to 0) and there is no charge contained in R = 0. In this case, 
( pUD can be expressed as^ 

As a matter of fact, the requirement of absence of singularities at R = is not 
necessary and, as is clear from the previous reasonings, one only needs the fulfillment 
of 

\imR 2 (m, R -R X (3, R ) = 0. (33) 

Thus, there is a stronger version of proposition |3.1| where the demand of absence of 
singularities is replaced by the milder condition (^), in which case there may be a 
curvature singularity. As a particular case of this, suppose that we can carry out 
Taylor's expansions of m and (3 in a neighborhood of R = 0: 

m(u,R) = mo(u) + mi(u) R + rri2(u )R 2 + 0(R 3 ), 
P(u,R) = p (u)+p 1 (u)R + p 2 (u)R 2 + O(R 3 ), 

then 

\q{u^ ^ ^ mi («)+4m ( M )A( M ) (34) 
where, as before, the righthand side is real due to fl2"T|). 



4 The Vaidya-Reissner-Nordstrom spacetime 

As stated in the introduction, the spherically symmetric space-time V will be taken 
as the Vaidya-Reissner-Nordstrom ( V-RN) solution [ E2\ , whose line-element is given 
by (D with 2m{u, R) = 2M(u) - Q 2 /R and (3(u, R) = 0, that is 

ds 2 = _L_ 2M^) + <g\ du2 + 2£ du dR + r2 dn2 ^ (35) 
y R R J 

where M(u) depends only on u and Q is a constant. Among the most remarkable 
known particular cases of ([$5|) we can cite: if M{u) is constant and Q ^ 0, then it 
coincides with Reissner-Nordstrom's (RN's) solution, including its particular cases 
of Schwarzschild's solution (Q = 0) and flat Minkowski's spacetime (M = Q = 0), 
2 (|27|) guarantees that the square root is real. 



4 THE VAIDYA-REISSNER-NORDSTROM SPACETIME 



13 



all of them very well-known and studied (see, e.g., |22|); if M{u) is not constant 
and Q = 0, then (R^) reduces to Vaidya's radiating solution |52|, whose structure 



and possible extensions have been treated in |33|, [T7| , including the possibility of 
producing naked singularities (see, e.g., @.) A general study of this solution was 
also carried out in |nj • 

Whenever Q ^ 0, the 2-form 

F = %- u A n . (36) 

satisfies Maxwell's equations in the metric (|35|) and thus Q is interpreted as a 
charge and (|36|) as an electromagnetic field present in the spacetime. The energy- 
momentum tensor T associated with fl3~5|) is given by 

T = -i_»»l + E (37) 
4:TiR 2 du 

where the electromagnetic part is given by ( |l"3D with $ = ip = Q/R 2 , clearly describ- 
ing a pure spherically symmetric non-null electrostatic field of total charge Q, and 
1 = — du is the future-pointing radial null vector tangent to the radial null geodesies 
with u =constant. Thus, the non-electromagnetic part of the energy-momentum 
tensor describes a kind of incoherent radially outgoing (respectively ingoing) radia- 
tion filling the exterior region of the spherical body whose total mass-energy, given 
by M(u), decreases with the retarded time u for e — — 1 (resp. increases if e — +1) 
as a result of the isotropic radial emission (resp. absorbtion). By analogy with (|14D , 
one might think that this kind of radial radiation can be in fact electromagnetic 
radiation. However, we have proved in section [3] that this is not possible, since 
there can be no null electromagnetic field without sources satisfying the Maxwell 
equations (|21~l - [24l ) in this space-time whose stress-energy tensor matches the form of 
T — E in (|37|), as should be expected 0. 

Observe that the energy conditions for the radiative part of the energy- momentum 
tensor hold if 

dM(u) , , 

e—^ > 0, 38 
du 

and that the V-RN spacetime is an example without the properties considered at the 
end of the previous section, because it has a strong curvature singularity at R — 0, 
and the milder condition fl33|) is also clearly violated whenever Q ^ 0. The global 



structure of this space-time if Q ^ was analyzed in [17]. There appear three cases 
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depending on whether or not there exists u\ such that M 2 {u\) = Q 2 and, if not, 
on whether M 2 (u) is greater than Q 2 or not. We refer the reader to ||T7| , |5T|| for an 
account of the particular features in each case and for the corresponding Penrose 
diagrams of the V-RN spacetimes. 



5 The gravitational-electromagnetic matching 

In order to match a general spherically symmetric space-time V with the Vaidya- 
Reissner-Nordstrom space-time V across a time-like hypersurface S, we are going 
to impose the gravitational and electromagnetic junction conditions. Specifically, we 
demand that the Einstein equations hold, in the distributional sense, in the whole 
space-time in such a way that no infinite jumps are allowed for the stress-energy 
tensor (thus, the possibility of a thin shell of matter on E is not considered). On the 
other hand, the Maxwell equations are also assumed to be valid, in a distributional 
sense, in the whole spacetime but allowing in this case for the existence of a charge 
surface density so that the electromagnetic 4-current may have an infinite jump at 
E. 

The practical technique to perform the matching is to consider V and V as 
divided, each one, by the candidate matching hypersurfaces a (respectively, a) into 
two parts -say 1 and 2 for V and 1 and 2 for V-. If a and a are diffeomorphic, one 
can identify corresponding points in them and then try to perform the gluing which, 
in principle, can be done in four different ways, namely, 1 — 2; 1 — I; 2 — 2; 2 — I. 
We can select one of them by choosing the relative sign of the normal vectors to the 
matching hypersurfaces a, a, see [|18j for details. Notice that, as these normal vectors 



are spacelike then we can always choose the tetrad such that n and n coincide, on 

a and a respectively, with them. Thus, from now on the unit normal vectors will 

also be denoted by n and n. 

Now, once a and o have been identified in the new glued spacetime there is no 

need to distinguish them and they will be simply termed as E if there is no confusion. 

If E is timelike and preserves the spherical symmetry then it can be described by 

intrinsic coordinates {£, y?} where £ is a timelike coordinate and where {«(£)> -R(0> 

9 = = ip} and {«(£)> -R(0> = 4> = are the parametric representations 

3 We will set an overbar on all variables, parameters and functions of the space-time V to 
distinguish them from the corresponding objects of the general space-time V. 
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of cr in V and of a in V, respectively. Before going any further, and as pointed out 
in |36], one has to specify how the tangent planes at every point p e o and at 
its corresponding point p G a must be identified in order to construct a well-defined 
geometry in the whole glued space-time. To do that, we proceed as follows: firstly 
we identify the vector fields associated to the angular variables, 8/89 with 8/88, 
and 8 /dip with 8 /dip, at E; secondly we consider the two vector fields defined on o 
and a respectively by 



(8R{£) 8 du{£) 8 



\ 8^ 8R 8i 8u t 



and 



(dR(£) 8 + du(£) 8 



(39) 



V 8£ 8R 8£ 8u / 

and identify them since they both represent <9/<9£ at E; thirdly, take the unit normal 
vectors n (resp., n) for a (resp., a) defined, except for a sign, by being orthogonal 
to the three previous vectors at E, and choose these signs in such a way that every 
curve crossing E through a point p = p must have a unique well-defined tangent 
vector there |Tj|; and fourthly, by doing so, the relative sign e n (e„ = 1) of the 
normal unit vectors has been fixed so that we can identify them thereby achieving 
a complete identification of the two tangent planes. As a by-product, this process 
determines an orientation and an arrow of time for the glued space-time if these 
concepts are well defined or fixed in some way in V and V. 

At this stage, we impose the Darmois gravitational junction conditions |13|||18||, 
which are the best suited for our purposes, requiring that the first and second fun- 
damental forms of E be identical when computed from either V or V . After the 
appropriate calculations, these matching conditions can be written in two (equiva- 
lent) complete sets, depending on the sign e n , as follows 

Case ee = e„ 



2 ^y 



Case ee 



H 



m, R 



2R 



R 

Eiie 2f3 
rn 

z 2/3 u + /3, R R 



R, 

EU, 

M — 



2R' 



£ 4R* U 



(40) 
(41) 

(42) 
(43) 



EX& u 



R 
2R 



R 



(44) 
(45) 
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m, R \ m, u e 2/3 



R 



m = M-^t, (46) 
e^u+(^-xi3,n)R I e^xfi. (47) 



Here = means that we have to compute both sides of the equality at S, and the 
overdots stand for derivatives with respect to £. 

In practice, we can understand the equivalence of the two sets of matching con- 
ditions by noting that they are suitable to describe the prolongation, when crossing 
S, of a given radial null geodesic with u =constant by either a radial null geodesic 
with u =constant (if ee = e n ), or not (if ee = — e n ), see |51] for more details. 



The above takes care of the gravitational matching, but we still require an ap- 
propriate junction of the electromagnetic fields. To that end, we present now the 
general matching conditions for the electromagnetic field (for arbitrary space-times) 
and then we will first particularize them to the spherically symmetric case and in 
a latter step to the case in which one of the matching space-times (V) is the V-RN 
solution. Let F and F be the respective electromagnetic 2-forms at V and V, so that 
they satisfy the Maxwell equations 

dF = 5F = 4ttJ (48) 

dF = 5F = 4ttJ (49) 

on V and V, respectively. The whole electromagnetic field is defined to be 

.F=F(1-0 E )+F0 S (50) 

where the Heaviside theta function 6*s is defined by #s|v = an d #s|v = 1- Thus, 
it is immediate to get 

dT = (FAn-FAn)5 E (51) 
5T = 4tt[J(1 - S ) + J# s ] + [z(n)F - i(n)F] 5 S (52) 

where A is the exterior product and 5 E is the normalized Dirac delta with support 



on S, see e.g. [12[ |36|], defined by d9% = n5 E . It follows from (|5T| , |52|) that the 
fulfillment of the Maxwell equations for T in a distributional sense provides the 
sought-after electromagnetic junction conditions: 

F A n = F A n (53) 
J = J(l-e s )+Jfe + O s , (54) 
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where J is the distributional 4-current and 

/cI^[i(n)F-z(ii)F] (55) 

is the surface current density 4- vector (obviously defined only on E). 

If we do not wish to allow the total 4-current to have infinite jumps at E, then 
from (^4j) and (|55"D we must impose 

i(n)F = i(n)F (56) 

so that, taking into account fl5H|) we derive F = F in this case. 



Let us now consider the particular case we are interested in: spherical symmetry. 
From conditions (|53|) one immediately infers that a non-null electromagnetic field 
(resp., a null one) can only be matched with another non-null field (resp., null). 
But as we proved before, null electromagnetic fields are not compatible with spher- 
ically symmetric energy-momentum tensors, so that only the non-null case must be 
considered for our purposes^. Hence, the electromagnetic fields take the form (|TT|) 

F = ipuAn + (f)UJgA uj v and F = ?p u A n + <j) u>g A u> v 

in V and V, respectively, so that the matching conditions ( |5*3|J>3| ) imply 

= (57) 
£ = i-$-V)u. (58) 

Let us finally particularize to the matching with the V-RN solution. Condition 
( pTD applied to the electromagnetic fields ( |Tl~D and (|36| ) becomes 

= (59) 

which, together with fllTD , leads to vanishing of on the entire V, and not just on 
E: 

= 0. 

This result means that the radially moving observers cannot measure any magnetic 
fields on V. Similarly, taking into account that the "charge function" on V is given 



4 



If one does not demand the invariance of the null field energy-momentum tensor and computes 



the junction conditions (p^,p3) for the null-null case one easily gets JC = and P = P, see 151 
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by q(u,R) = ipR 2 [see (P0|)1, we can define the iota/ charge Q = ?/>-R 2 . Using now 
fl36|), we can express /C of fl58]) as 



so that, using 



we are led to define the total surface charge in E by 



This is a very natural result and implies that the sum of the total charge Q con- 
tained in V plus the total surface charge must be constant (= Q). Of course, 
the presence of infinite jumps on the electromagnetic 4-current is a mathematical 
idealization that can, in many cases, simplify the complexity of describing some 
models with a large concentration of charge at E. However, the absence of infinite 
jumps in the 4-current is more realistic physically and, in fact, it is in agreement 
with the experimental observations J2"8| . Therefore, if in particular we do not allow 
for infinite jumps in the 4-current or, in other words, if there is not a surface charge 
we obtain 

(60) 



Q = ^R 2 = Q. 



Then the total charge Q must also be constant. 

In the rest of this section, we are going to focus on this physically more realistic 
case when (|6T)D is fulfilled. Then, the complete set of matching equations fl4"0|)-(|4*3|). 
( |59"|) and flBTf ) (or alternatively ([44D -(f47|), and (^3)) are the necessary and suffi- 
cient conditions for the matching of a general spherically symmetric metric with the 
V-RN solution through a general spherically symmetric timelike hypersurface, and 
they provide relations between relevant quantities at both sides of E. Nonetheless, 
it is remarkable that a linear combination of (f4"0"D, (|B"0"D and (E3) leads to 



2m\ 



R J y,R 2R ' 4R : > 
or alternatively, for the other case with ee 



m, R Q 2 



x(xP, 



2R 



+ e 



m m e 
R 



+ 



4# 3 



u+\ 



u-f3, R e- 2 ?R^0, 
n , one similarly finds 

Q 2 



(61) 



m, R 
R 



R 



2R 3 



e" 2 ^ = 0. (62) 



Clearly, these relations involve quantities of the space-time V only, but not of V. 
This means that equation (|6l|) (or (|62j)) is a necessary condition that the space-time 
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V, by itself, must fulfill on S in order to be matchable to a V-RN solution. The 
physical meaning of (|6lD (or fl62|)) generalizes the results for standard models in 
which normal pressures must vanish on E when matching to a vacuum, as will be 
explained in section in detail. 

Given the previous remarks, one can now derive much more information from 
the matching equations. To extract this information, one has to follow paths which 
depend on the known data of the problem under consideration. A physically very 
interesting case arises when m(u,R) and j3(u,R) are given, so that the spacetime 
V is completely known, and we want to ascertain whether or not V is matchable 
to a V-RN solution, if so where, and finally to which particular V-RN metric, that 
is, for which particular M(u) and Q. In order to solve this problem, first of all 
we treat Q(— Q) either as known or as a parameter on which the final results will 
depend. This is a key point. Then, we proceed as follows: for the ee = e n case 
(resp., ee = —e n ), the equation (|6l|) (resp., (|62|) ), can be considered as an ordinary 
differential equation for R(u), with the form dR/alu = F(u, R; Q), so that if m(u, R) 
and f3(u, R) are such that F(u, R; Q) satisfies Lipschitz's conditions we can find the 
solution R(u;ci,Q), where c\ is an integration constant. If R(u;ci,Q) defines a 
time-like £ on V, then using (fEq) we can determine M = M(u;ci,Q). Now, by 
integrating (|4l|) for the first case -or (f|^) for the second- we get u(u; c\,Q), except 
for a new additive constant C2,Q while from ( fl0|) (or, respectively, from fl44|)) we 
get R(u; ci, Q). These two functions define the hypersurface £ in the space-time V. 
Finally, due to the fact that M(u) is a function of only u, by combining M(u; Ci, Q) 
with u(u; Ci, Q) we get M(u; ci, Q) and the problem is completely solved. 

The solutions R(u; ci, Q) give a two-parameter family of matching hypersurfaces 

through which the space-time V is matchable with a V-RN solution. The different 

values for c\ (and of course Q) will simply lead to the different particular V-RN 

space-times which match at the various S's, providing explicitly M(u) (and Q). In 

order to interpret the physical differences between these S's, consider their intersec- 

tion with a given null hypersurface u = uq =constant. The values of R = R and 

M at these intersections are denoted by R = R(u ; c±, Q) and M = M(uq; Ci,Q), 

so that they depend on the value of Ci, that is, on the particular E. Eliminating c\ 

from M and Ro we obtain M = M (uo, Rq,Q) which gives, at the given uq, the 

5 The specific value of C2 is irrelevant since it only provides an origin for the times u and u, and 
therefore we choose, for the sake of simplicity, ci = 0. 
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value of Mq as a function of Rq. Use of fl42|) [or (f46|)1 leads to the explicit expression 

Q 2 



M (w , i? , Q) = m(u 0: R ) + 



2R 



o 



Differentiating here with respect to Rq and taking into account the energy condition 



( |27D with Q = Q we get 



2 



Ml 



Finally, if the hypersurface at w is not at a region with trapped 2-spheres, that is 
to say, x( u 0y Ro) > 0, then from the energy condition @ we finally infer 

dM 

> 0. 

dRo 

In other words, given a space-time V take the 2-spheres (which belong to a matching 
hypersurface S) defined by £ fl {w = w } H {x > 0}. These will have an area 
proportional to R^ and contain a total mass M . The previous result tells us that 
smaller such 2-spheres will contain non-bigger masses. 



6 Implications on the matter content of V 

In the previous sections we have not interpreted the matter contents of the space- 
time V apart from assuming that it may contain an electromagnetic field. Our 
purpose in this section is to investigate some of the different interpretations for the 
non-electromagnetic matter content and to elucidate the restrictions imposed by the 
matching conditions, found in the previous section, on the electromagnetic field as 
well as on the non-electromagnetic part of the energy-momentum tensor. 

Let us start with the electromagnetic field in V, given by ([LJ]) with ([18]) and 

ij = Q/R 2 , = 0. 

Taking the derivative of the first of these equations along S we get, on using (|18|) 

i{u)d{R 2 i)) = J = ^-^i(n)d(R 2 ij)u. (63) 

This tells us that J is a time-like vector on S, furthermore comoving with £ at 
S. Outside S, though, the form of J is unrestricted, and in general it has two 
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components, one along the direction of the time-like u, which is the convection 
current with respect to u (assuming that u is chosen to be the 4-velocity of the 
charged fluid), and another in the direction of n, representing the corresponding 
conduction current with respect to u, which necessarily vanishes at S. Therefore, 
any function ift that satisfies flBTf ) provides us with a regular electromagnetic field 
valid for V.Q 

Now we can ask ourselves the physically relevant question of whether the electro- 
magnetic field can be chosen only with convention, or only with conduction, current. 
This is clearly equivalent, according to fll8|) , to finding out the causal character of 
the one- form d(R 2 ip), because the sign of J • J is opposite to that of d(R 2 ip) ■ d(R 2 ip)- 
Thus J is time-like (resp. space-like) and thus can be aligned with u (resp. n) if 
d(R 2 ip) is spacelike (resp. timelike). Observe that in the second case J = neces- 
sarily. 

Finally, the energy-momentum tensor E associated with the non-null electro- 
magnetic field in V is now (|T3|) with = 0: 



ij 2 

E = — (u<g)u - n<g>n + u) e ®u) e + u;„<g)u;A (64) 

57T 

Let us now pass to the rest of the matter content. Hitherto, we only know 
that its energy-momentum tensor, P = T — E, takes a form of type (|8|). As is 
known, see e.g. [f48j, the energy-momentum tensors for different physical types of 
matter distributions can, in fact, have precisely the same components. Therefore P 
is opened to different interpretations, depending on the physical processes we want 
to describe. 

One of the simplest interpretations that we can put forward, containing the 
minimum degrees of freedom, is the case of a charged non-perfect fluid embedded in 
a radially directed null radiation. The 4-velocity is taken to be u so that the fluid is 
comoving with E by construction. The total energy-momentum tensor for this case 
is then written as 

T = (fi + p) u®u + p g + tt 2 £®£ + n + B, (65) 

where \i is the energy density of the fluid, p its isotropic pressure, II its anisotropic 
pressure tensor (which is traceless and orthogonal to u), Q 2 is the null radiation 
6 If wc allowed for infinite jumps on J the problem would be easier since there would be no 
constraints on ip, not even on the matching hypersurface, so that any function tp would provide a 
regular electromagnetic field valid for V. 
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energy density, £ = u + s n (s 2 = 1) is a radially directed null vector field, and E is 
given in (|64]). 

Obviously, the form (^) of the energy-momentum tensor is not unique and 
another physically relevant possible interpretation, also having the minimum degrees 
of freedom, IS cLS cL charged non-perfect fluid with heat conduction. In this case we 
write 

T = (/i+j5)u<g>u + pg + h<g>u + u<g>h + n + E, (66) 
where h oc n is the heat flux vector and jl, p and II have the usual interpretation. 



Of course, (p5,66) are not the only possible cases and many other physically realistic 



decompositions can be written according to the particular situation. 



In any case, and following the philosophy proposed in fllBI , for every given 
m(u, R), f3(u, R) and ip(u,R) such that (|60|) is satisfied and a timelike £ exists, 
the velocity vector field u is fixed only on X and, provided that it is time-like and 
orthogonal to the 2-spheres {u, R} =constants, we have freedom to choose it in the 
rest of V. Once u is chosen we can compute, by means of ([Tj|), the 4-current J, and 
by using Einstein's equations, /i, p, Q 2 and II for (|65D, or fl, p, h and II for (|66D , or 
the corresponding quantities for any other desired interpretation of P. 

It is known that some relations between the energy-momentum tensor at both 



sides of S follow from the matching of two space-times PH B61. They read 



T^n" = T F fY, (67) 



- [IP 



n^n v = T^n v , (6£ 



where r M (equivalently r M ) is any vector tangent to E. Therefore, for the first 
interpretation proposed in (|65D, the non-trivial relations deriving from ( |67| , |68D are 



I = 1, (69) 



where we have used ( |S^) and d37|) . Combining the last two equations and using (160| ) 

we get 

P + n M x^ = o, (72) 

which does not involve any quantity from the space-time V. 
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The equalities (|60D, (|70|) and (|72|) are the main physical equations which relate 
quantities at both sides of S. The equation (fTCf ) expresses the fact that the radiated 
energy density has to be continuous through the matching hypersurface. On the 
other hand, it is easy to check that ( |72|) is strictly equivalent to the matching con- 



dition (plD (or to (p2|)), see Its physical interpretation is clear, (^) says that 



the total normal pressure has to vanish on the matching hypersurface. All in all, we 
can reformulate the results from this and the previous sections in the following way 



Theorem 6.1 Every spherically symmetric metric can be locally matched to a Vaidya- 
Reissner-Nordstrom solution across any time-like hypersurface with the properties: 
1 ) the total normal pressure of the fluid in the decomposition ME) vanishes on it; 
and 2) the total charge enclosed in it is constant. 

This is a very satisfactory result, which generalizes previous results for the uncharged 
case, in particular those obtained by Misner & Sharp [^] and Bel & Hamoui for 
the case with a Schwarzschild exterior and absence of radiation; and those of Fayos 



et al. [16 1 for the Vaidya's radiative solution. It is also a generalization of previous 



results obtained by Oliveira & Santos j|2| for charged shear-free fluids, although this 
case is better adapted to the second interpretation given in (|66|). We can actually 
reformulate the previous theorem in this case as follows: 

Theorem 6.2 Every spherically symmetric metric can be locally matched to a Vaidya- 
Reissner-N ordstrom solution across any time-like hypersurface whose total normal 
pressure in the decomposition j^Udj ) equals the heat flux on it, and such that the total 
charge enclosed inside it is constant. 

6.1 Brief note on the existence of a surface charge on £ 

If we allowed for the existence of infinite jumps in the electromagnetic 4-current the 
situation would be slightly different, and for example condition flT2|) would not be 
valid. Instead, we would have, from (|7(]) and (|7T|): 



Q 2 e Q 2 

" 8-/." 8-/."' 



p + n^n^n" 



This can be easily understood by noting that the normal tensions exerted by the 
electromagnetic fields on S, given by —Q 2 /(8nR 4 ') and — Q 2 / (87iR 4 ), do not com- 
pensate reciprocally. A non-zero "normal pressure exerted by the fluid" is thus 
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needed in order to satisfy the equality of total normal pressures as required by fl68D . 
Given this, we can generalize theorems |6.1| and |6.2| as 

Theorem 6.3 If infinite jumps in the electromagnetic ^-current are allowed, ev- 
ery spherically symmetric metric can be locally matched to the Vaidya-Reissner- 
Nordstrom solution across any time-like hypersurface £ such that 

AttR a (p + U^n") - Q 2 = constant < . 



Theorem 6.4 If we allow for infinite jumps in the electromagnetic 4-current, every 
spherically symmetric metric can be locally matched to a Vaidya-Reissner-Nordstrom 
solution across any time-like hypersurface satisfying 

4vr J R 4 (p + U IM/ n li n u - h a n a ) - Q 2 = constant < . 

Clearly these theorems are less restrictive than the corresponding theorems in the 
absence of surface charge, since the time-like hyper surf aces satisfying the require- 
ments in theorem |6J] (or |6.2| ) will also satisfy the corresponding requirements in 
theorem |673] (or, respectively, |6.4|) , but not the other way round. 



7 Other important physical consequences 

In this section we are going to prove that the energy and matching conditions imply 
several important physical consequences, namely, (i) inequalities valid on the match- 
ing hypersurface, one of them involving only quantities of V; (ii) the non-spacelike 
character of the matching hypersurface under general conditions; and (iii) several 
limits on the total charge of the model. 

To that end, we are going to exploit the energy conditions on P, given by (|3]), 
(|5|), and (p7|-|2"9"l) supplemented with the fact that = and, therefore, $ 2 = ip 2 . 
This group of conditions constrain the functions m(u, R), /3(u, R) and ip{u, R) that 
we can choose for V and, when specialized to S, give physical restrictions valid on 
the matching hypersurface. For instance, the matching conditions for the ee = e n 
case[] are the relations (|40|)-(42) and (0). We can isolate m t R and m jU from them 



7 Obviously the ee = —e n case can be treated similarly. 
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in terms of magnitudes from V 

s 



m, R = 2RUx-eR!)^ R + ~\ (73) 



em, u I ee 2 "(M -2^ iJi (x-e^) J R') (74) 
where i2' = .R/w. Combining these with ( p7| ) and (|5]) we arrive at 

< RP, R {X-^R') (75) 
< eM fi < l -Rf3 tR {x ~ 2eR') 2 (76) 

where use have been made of the energy condition (^). On the other hand, if we 
differentiate ( f£2"D on £ and replace Rjix with the value arising from (^T|) we find 

£m, u (3, R > —(Q 2 - 2R 2 m, R )(Q 2 - 2R 2 m, R + A X (3, R R 3 ). (77) 

where again the fulfillment of (^) has been used. Unlike ([F^) and (|76|), this in- 
equality (|77|) only relates quantities from V on a, and it will be important in the 
construction of models. On the other hand, inequalities ([T§), ( [76] ) and ([H]) use up 
all the reciprocal implications of the dominant energy condition at both sides of the 
matching hypersurface. 

More interestingly, we can derive under which circumstances the matching hyper- 
surface will be actually non-spacelikeQ. The condition for this is u(2eR — e 2l3 x'u) < 0. 
We can consider two possibilities: 

1. it = 0, so that £ is null. Then we get: 

• Case ee = e n , from fl61~|): 

P,R=0 



Case ee = — e n , from (p2|) : 

s 



R 2R 3 



2. u > 0, so that 



R e e 2/3 
u 2 



Isolating Rjix from (|6T| ) [or from (^) for the second case] and replacing it into 

we get 



8 Here we will use that (61) and ( p2| ) are valid no matter what the character of the matching 
hypersurface is, as was shown in M. 
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Case ee = e r 



H WO) =► ^-^-^>0 (79) 



CclSG — 6^ '. 



I if ~^'-R~2^3 ^ °) ^ ^ >0 (80) 

These relations are satisfied if ( p7|) and (|5]) hold. Furthermore, > and 
X 2 ft,R + e (2^, u e~ 2/3 )/R > 0, wnt/j £/ie s£hc£ inequality, are part of the energy con- 
ditions for energy- momentum tensors of type I (see |22| for definitions). Therefore, 
we have 

Proposition 7.1 7/ the energy-momentum tensor in V is of "type I" and its non- 
electromagnetic part satisfies the dominant energy conditions, and if no infinite 
jumps are allowed for the electromagnetic ^-current, then any matching hypersurface 
£ with a Vaidya-Reissner-Nordstrom solution is necessarily non-spacelike. 



More specifically, assuming that ( f2~7[ ) holds, it is easy to see that the matching hy- 



persurface will generically be timelike, and it may be null only in the very particular 
case that m, R /R - xP,R -Q 2 /(2R 3 ) = holds. 

On the other hand, type II energy-momentum tensors are usually associated 
with zero rest mass fields (recall that, for instance, the Vaidya and V-RN solutions 
have a type II energy-momentum tensor). We can consider the question of whether 
or not a proposition equivalent to |7.1| (relating the energy conditions for V with 
the non-spacelike nature of S) can be found when the energy-momentum tensor is 
type II. It is easy to see that such proposition cannot exist in general, since the 
matching conditions in the previous section show that, in the particular case of 
two matchable V-RN space-times, the nature of the matching hypersurface can be 
arbitrary independently of any other conditions. 

Let us consider finally the important question of whether there are any upper 



limits on the total charge of the model. From proposition |2.2| evaluated on E and 
using (|42]) [or (f46|)1 we immediately obtain 



Proposition 7.2 lfm(u,0) > and the dominant energy conditions hold, then 
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This is a key result in our treatment, for it can be interpreted as a constraint on the 
charge enclosed in S. Notice that this constraint exists even if there are curvature 
singularities at R = provided that m(u, 0) > 0. Ponce de Leon [J46] (see also 0) 
obtained a similar result for the much more restricted case of static, singularity- 
free space-times with a charged perfect fluid matched to a pure Reissner-Nordstrom 
exterior. A general result for the general class of static singularity-free metrics can 
also be deduced from the limits found in [53], and a more recent work is presented in 
[ P^ ] . The fact of considering only static space-times is, possibly, the most important 
restriction in the above-mentioned papers, since then the results can only be applied 
to "limiting configurations" for charged static spheres. Only our general approach 
allows for an interpretation of ( |STD in terms of the dynamics of the problem and 
allows for a general analysis of the possible existence of naked singularities (see e.g. 
|S| jf29] for definitions and relevance) and their creation. 

With regard to this, observe that (|8~l]) can be also viewed in other interesting 



ways when expressed as 

Rz(u)> Q 



2M(u)' 

This inequality implies that any space-time V providing a physically realistic model 
for the interior of a collapsing star either has its limit surface at values of Rz(u) > 
Q 2 / (2M(u)), or otherwise a necessarily timelike singularity^ (therefore locally naked) 
must develop at R — 0. The question on whether this radius can effectively be 
reached lies beyond the purposes of this article, but you can consult |50|| and refer- 
ences therein for an interesting discussion at this respect. 

There are also constraints on the charge Q of V due to the the fulfillment of 
DEC. One of them is given by (]77|). Apart from this, inequalities ( |2"71 ) and (|2"8"| ) 
define two new possible maximum values for Q: the first one is just the appropriate 
specialization of ( |30]) 

Q 2 < 2R 2 (m, R -Rxf3, R ), (82) 

while the second one reads 

Q 2 < l -R\\ + P2). (83) 

The righthand sides of these two inequalities vanish at R — if there is no curvature 
singularity there, as is clear from ([|). Thus, if S arrived at the vicinity of R = 
with non-zero charge, the previous conditions would be violated from a certain value 
9 The singularity must be timelike because m(u, 0) < 0. 
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of R on (decreasingly). This result is similar to that in proposition \l.2\ S can only 
get arbitrarily close to ill = if either a curvature singularity already exists, or 
develops, there. 



8 Conclusions 

The aim of this work is to provide a theoretical framework for the construction 
of global models describing stars and voids in General Relativity. Our only as- 
sumption is, apart from the spherical symmetry of the spacetime, that the exterior 
of the star (respectively the gravitational field within the void) is represented by a 
Vaidya-Reissner-Nordstrom solution properly matched to an interior (resp. exterior) 
spacetime with the only restriction that the dominant energy conditions are fulfilled 
everywhere. Thus, we are treating the general case in which the stars/voids may be 
charged and radiating. 

A summary of our main results is the following. To start with, in proposition 



272] we saw that the DEC imply the non-negativity of the mass function everywhere 
if m(u, 0) > 0. This result has important consequences in the rest of the paper, 
as for instance the restriction which imposes on any spherically symmetric electro- 
magnetic fields, see section |3]. More important consequences are the restrictions on 
the presence and distribution of electric charge, first in a neighbourhood of R = 



(section 0) and then on the whole spacetime (formulae (]g2"|-[g5D and proposition |7.2|) 



As a specially important consequence, we have shown in full generality that the 
total charge of any physically acceptable spherically symmetric object has un upper 
bound related to its size and mass given by Here, by physically acceptable we 
mean that the object satisfies m(u, 0) > which in particular includes all bodies 
which are regular everywhere. This is a result genuine to General Relativity for, as is 
well known, Classical Electrodynamics does not impose any limitation on the charge 
that an object can possess. Furthermore, our result is completely general and does 
not depend on the assumptions of staticity used in its previous partial and particular 
versions found in EE], |5], and references therein, or in the more recent ||34j| . In 



this last work an inequality relating the mass, the radius and the charge of general 
static matter distributions has been found by using the gravitational field equations 
(in a similar way as Buchdahl obtained the classical maximum mass-radius ratio 
for uncharged stable stars). Their result shows that a charged (Q < M) stable star 
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must obey the inequality (fH]). However, as is clear, our treatment provides the 
limit ([81]) in dynamical situations, including evolving voids or the important case of 
collapsing, or rebounding, stars. 

The appearance of a maximum value for the charge can be alternatively read as 
the existence of a minimum value for R%(u). This result can be interpreted, for the 
case of a charged star -no matter how small its charge is- that obeys the DEC as 
stating that, if the star ever reaches values of R lower than the minimum radius, 
then a necessarily timelike singularity, ergo locally naked, has to develop at R = 0. 

All our results have been derived by properly matching the interior and exterior 
spacetimes, and here we have taken into account the pure electromagnetic matching 
conditions as well as the gravitational ones, which were fully derived in section [5] 
both in the case with no infinite jumps for the total electromagnetic 4-current - 
which agrees with experimental results — or in the case with them — this is the usual 
approach in the literature but, in our opinion, it is unrealistic from a physical point 
of view, except as an approximation to the real world — . We have also proved that, 
in generic situations, there is a 2-parameter family of matching hypersurfaces which 
were later shown to be non-spacelike in Proposition |7!T]. The general result that 
the vast majority of spherically symmetric spacetimes can be matched to a V-RN 
solution (or one of its particular cases) has a clear physical interpretation given by 



Theorems rO and |6.2|. These results generalize the ones obtained in |T6l 39 



for several particular cases arising if there is no charge or radiation. However, let us 
remark that even if the total charge of the modelized object is zero and, consequently, 
the metric for the space-time V is the Vaidya metric (or its specializations), still a 
radial distribution of charge is allowed in V, even a time-dependent one, as long as 
the total charge is zero. This fact allows for a possible description of the charge 
redistribution phenomena which could take place in the interior of a star. This is 
the case of, for instance, the hybrid stars, in which, despite global charge neutrality, 
the local non-neutrality is energetically favoured (see pOj ) . 

In our next article (Spherically symmetric models for charged radiating stars and 
voids II: Practical approach) we analyze in depth some of the theoretical formula- 
tions proposed in this work and examine some models for stars and voids under the 
scope of our approach. 
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